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Distributed source coding (DSC) addresses the compression of correlated sources without communi- 
cation links among them. This paper is concerned with the Wyner-Ziv problem: coding of an information 
source with side information available only at the decoder in the form of a noisy version of the source. 
Both the theoretical analysis and code design are addressed in the framework of multi-dimensional nested 



C/3 ' lattice coding (NLC). For theoretical analysis, accurate computation of the rate-distortion function is 

given under the high-resolution assumption, and a new upper bound using the derivative of the theta series 
is derived. For practical code design, several techniques with low complexity are proposed. Compared to 
[> — ' the existing Slepian-Wolf coded nested quantization (SWC-NQ) for Wyner-Ziv coding based on one or 



two-dimensional lattices, our proposed multi-dimensional NLC can offer better performance at arguably 
lower complexity, since it does not require the second stage of Slepian-Wolf coding. 

Index Terms 

Distributed source coding, multi-dimensional lattices, nested lattices, rate-distortion function, Wyner- 
Ziv coding. 

I. Introduction 

In wireless sensor networks, energy is a major constraint to individual node performance. 
Compressing data prior to transmission and minimizing inter-node communications can effec- 
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tively reduce energy consumption. The problem is well-known as multi-terminal source coding 
or distributed source coding (DSC) [[U, where one or more of the sensor nodes (or terminals, 
information sources) compress data separately (i.e. without communication with each other) 
before transmission. The joint decoder at the central node recovers the transmitted data either 
losslessly or up to a prescribed distortion. For the lossless case, Slepian and Wolf [|2| proposed 
the theoretical framework. A common scenario for the lossy case is the Wyner-Ziv problem, 
which considers encoding a source where the decoder has access to some side information [3J. 

Code designs for Wyner-Ziv coding have emerged in recent years. A structured binning scheme 
using (random) nested lattice codes (NLC) was introduced Hi, which can asymptotically 
achieve the Wyner-Ziv limit [3] as the dimension of the lattice approaches infinity. The first 
practical code design based on this idea was done in [61], where preliminary analysis for particular 
nested lattices was presented. A general distortion analysis for n-dimensional lattices was given 
by Liu et al. [7J. In particular, a lower bound on the distortion was derived, which, for given n, 
exhibits an increasing gap to the Wyner-Ziv limit as the rate grows. For this reason, |I71 focused 
on 1- and 2-dimensional lattices and used a second stage of coding, namely, Slepian- Wolf coding, 
to further exploit the correlation between quantized data. This technique was termed as "Slepian- 
Wolf coded nested quantization" (SWC-NQ) in fT]. The Slepian- Wolf coding in [TJ relied on 
channel capacity-achieving codes such as low-density parity check (LDPC) codes, resulting in 
considerable implementation complexity. 

In this paper, we take a new approach by directly using multi-dimensional NLC, which is 
conceptually simpler. It was known that the distortion performance improves with the lattice 
dimension JH, fTl. Our work is inspired by this result, and we aim to develop less complicated 
practical codes for the Wyner-Ziv problem. The contributions of this paper are two-fold. On one 
hand, we complement the theoretical analysis of the rate-distortion function in [7J. Specifically, 
we develop a technique to compute the distortion, which is accurate under the high-resolution 
assumption, and a new upper bound based on the derivative of the theta series. It is worth 
mentioning that, in practice, an upper bound is a safer guideline of system design than a 
lower bound. On the other hand, we implement NLC based on a variety of techniques to 
obtain sublattices, and show that our implementation gains over the conventional one and two- 
dimensional schemes. These techniques include clean similar sublattices of dimensions greater 
than two [SJI, the random ensemble of nested lattices in [|9||, and scaling and rotation of a lattice 
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to obtain its sublattices. Recently, a similar idea using the scaling-rotating technique to obtain 
nested lattices was proposed independently in [[T0| for the synchronous multiple-access channel. 

The obtained distortion performance is close to the Wyner-Ziv limit, especially at rates < 
10 bits/sample. These simulation results are consistent with theoretical analysis which predicts 
improved performance as dimension increases, although the law of diminishing returns will apply 
in high dimensions. Meanwhile, the complexity of our proposed scheme is arguably lower than 
that of SWC-NQ, as we do not use turbo or LDPC codes El U^. 

11. Preliminaries 

Consider Wyner-Ziv coding in the two-source case. Let {{Xi,Yi)}^-^ be a sequence of inde- 
pendent and identically distributed (i.i.d.) drawings of a pair of correlated random variables X 
and Y, and let d{Xi, Xj) denote a single-letter distortion measure between the source Xi and its 
reconstructed version Xi at the decoder. Wyner-Ziv coding [[31 asks the question of how many 
bits are needed to encode source X under the constraint that the average distortion E[d{Xi, Xi)] 
is not greater than a given target distortion D, assuming the side information Y is available at 
the decoder but not at the encoder. 

We specifically consider the quadratic Gaussian case of the correlation model X = Y + Z, 
where Z is the Gaussian noise with distribution Z ~ A^(0,cr|), and Y is independent of Z. 
Mean-square error (MSB) is used for distortion measure. The quadratic Gaussian case is special 
since it has no rate loss, namely, the rate-distortion function is given by ATI 

Dwz{R) = 42-''' (1) 
which will be referred to as the 'Wyner-Ziv limit'. 

A. Lattices, nested lattices and similar sublattices 

For a set of n linearly independent basis vectors {mi, m2, . . . , m„}, an n-dimensional lattice 
A is composed of all integral combinations of the basis vectors: 

A = {1 = M • i : i G Z"} (2) 

where M = [mi, m2, . . . , m„] is the generator matrix and Z = {0, ±1, ±2, . . .} is the set of 
integers. For a vector x = [xi, X2, . . . , x^]^, the nearest-neighbor quantizer associated with A 
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is Qa(x) = argmiiiigA ||x — 1||. The basic Voronoi cell of A, defined by V = {x : Qa(x) = 
0}, specifies the nearest-neighbor decoding region. Important quantities for V include the cell 
volume V = dx, the second moment o"^ = ^ /y | |x| prfx and the normalized second moment 
G(A) = a^V-. The minimum of G(A) of all the n-dimensional lattices is denoted as 
From [fT3l . G„ > l/2ne,\/n and lim„^oo = l/27re. 

Let Ap he a fine lattice with a generator matrix M^. Similarly, let A^ be a coarse lattice 
with a generator matrix Mc. A pair of n-dimensional lattices {Af,Ac) is nested in the sense 
of Ac C A^, if Mc = Mi? ■ P, where P is an n x n integer matrix with determinant greater 
than one (in absolute value). We define the nesting ratio N = Vc/Vp, where Vp and Vc are the 
cell volumes of the fine and coarse lattice, respectively. 

There is one special case where Ac is geometrically similar to A^, which means that Ac can 
be obtained from A^? by applying a similarity transform lIH including a rotation, change of scale 
and possibly a reflection. Ac is strictly similar to Ap when reflection is not used. We also refer 
to Ac as a similar sublattice to Ap. 

B. Encoding and Decoding Scheme 

Throughout the paper, we follow the high-resolution assumption in [[7]|, which means Vp is 
sufficiently small such that the probability density function (pdf) of Z is approximately constant 
over each Voronoi cell of the fine lattice Ap. 

The encoding and decoding scheme is the same as in [7J, which is simplified from |[5l under 
the high-resolution assumption. The scheme is described as follows: 

• The encoder quantizes x to xg^^ = (5a^(x), compute the coset leader s = xg^^ — 
Qac(^Qa^)' ^iid transmits the index corresponding to s. 

• The decoder receives s and reconstructs x as x = s + QacIy ~ s). 

Under this formulation, the rate per dimension is given by i? = ^ log2(^). In our practical 
design, we also use minimum MSB estimation as described in [jVl. 
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C. Theta series and even unimodular lattices 

The development of this paper is heavily based on the theta series. Given a lattice A, the theta 
series [[T3l is defined as 



xeA 

where q = e^'^'^ (Q'(t) > 0). Using the change of variable jr — )• r (r > real), it can alternatively 
be expressed as 



Theta series of some standard lattices are given in [[T3l. In general, theta series are not easy to 
evaluate, and their derivatives are even more complex to calculate. Fortunately, there is a class 
of lattices for which the computation of such functions is tractable, the even unimodular lattices. 
A detailed description of such lattices can be found in [14]. 

A lattice A is unimodular [13] if A is integral, i.e., M'^M is an integer matrix, and | det A| = 1 
or equavlently A is equal to its dual A*, i.e., A = A*. If A is integral, then ||l|p is necessarily 
an integer for all 1 G A. Further, if ||lp is an even integer for all I G A, then A is called an even 
unimodular lattice. Many exceptional lattices such as the Gosset lattice or the Leech lattice 
A 24 are even unimodular. 

We introduce now the three Jacobi theta functions [[T3l . 
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Now, consider the Eisenstein series. 




(5) 



where C, (s) is the Riemann zeta function. 




fc=i 



We can relate (g) and the Jacobi theta functions through the relation. 
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Another fundamental series is the so-called modular discriminant defined as, 

A(?) = Y^(i^l(?)-^6(^)) (6) 
which is also related to the three Jacobi theta functions through the identity 

Remarkably, theta series of all even unimodular lattices can be expressed as polynomials in 
the two variables E,i and A: 

Proposition 1 ([fT5|). If A is an even unimodular lattice of dimension n then, 

1) n is a multiple of8,n = 24m + 8 A; with k = 0,1,2 and m being any non-negative integer; 

2) its theta series is given by 

m 

Oa (r) = (e-^^O + (e-'"0 A^' (e-^-^) , b, e Q. (7) 

i=i 

III. New Analysis of Rate-Distortion Function 

Knowing the suitability of lattices to this problem through theoretical analysis is essential for 
the code design. Our starting point is the rate-distortion function of nested-lattice Wyner-Ziv 
coding derived in |I71. It was shown there that the distortion D is comprised of a 'source- 
coding component' Ds and a 'channel-coding component' Dc- More precisely, the distortion 
per dimension for the Wyner-Ziv coding at high resolution is given by [17| 

Dn = mm {Ds + Dc} 

= min(G(A^)\/j + -Ez[||gAc(Z)ir] 
vcXJ I n 



= min |G(A^)v|2-2« + i \M'P{z e Vcih))] 
L iceAc J 



(8) 



where Ds and Dc are implicitly defined, and VcO-c) is the Voronoi cell associated with the 
lattice point Ic G A^. 

There is a standard way to handle Ds'- choosing a fine lattice that is good for quantization, 
i.e., a lattice with G{Ap) as small as possible. The second component Dc is not so easy to 
handle, though. This is the main problem we want to solve in this paper. 
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A lower bound on D„ was given in where the Voronoi region VcO-c) is replaced with 
the packing sphere S [fT3|| . It was pointed out in that this lower bound is asymptotically 
tight since the shape of VcO-c) will approach a sphere as n increases for lattices with the best 
rate-distortion performance. However, for finite n, the Voronoi regions may not be close to a 
sphere. We observe that for these lattices, the lower bound results in a non-negligible gap. 

A. Accurate Calculation Under High-Resolution Assumption 

We present an accurate calculation of the rate-distortion, based on the high-resolution assump- 
tion. 

Theorem 1 (Accurate Calculation Under High-Resolution Assumption). Let Ip G Ap be a fine 
lattice point and let 

be the pdf of a vector z whose elements Zi's are i.i.d. Gaussian with variance oz- Then, under 
the high-resolution assumption, the channel coding component Dq of the distortion function for 
nested-lattice Wyner-Ziv coding is given by 

E ii^^ii' E f^iw^F + icW). (10) 

Proof: Based on the high-resolution assumption [|71, the pdf /z(z) is uniform over each 
Voronoi cell of the fine lattice. Thus P(z G Vc(lc)) in ^ can be approximated by 

P(zeVc(lc))= E /z(I|1f + 1c||)V>. (11) 

Then the theorem follows. ■ 
To compute (flOl) . we need to enumeration all lattice points Ic G Ac, and those points If G Vc- 
The latter may be done by checking all points Ip within the covering radius of A^: li? G Vc if 
and only if \p is decoded to in Ac. 

In Fig. [H we give the rate-distortion performance of fine lattice and Eg. The sublattice 
is generated from the expand-rotate method described in the next Section. We also give the 
corresponding lower bound in ^ for comparison. From Fig. [T] we can see in dimension eight, 
the lower bound in is tight for the best quantization lattice Eg in this dimension, but has a 
gap to our theoretical calculation for the lattice. We find this is also true for other dimensions. 
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B. An Upper Bound 

Albeit accurate under the high-resolution assumption, the method in the preceding subsection 
unfortunately does not lend much insight into the selection criterion of a good coarse lattice. 
Let us reexamine the Dc component in ^ 

which should be minimized. Of course, this problem itself is not tractable due to the integration 
over the basic Voronoi region. 

To deal with this problem, the common approach is to choose a coarse lattice to maximize 
P(z G Vc), which is reduced to the standard approach to channel coding. However, this standard 
approach does not quite capture the essence of Wyner-Ziv coding, since the metric in (fT2l) is 
clearly different. Alternatively, one could derive either the upper or lower bound by using the 
covering or packing radius, the latter of which has already been done by Liu et al. LZ|. Again, 
in addition to the computational complexity, it is not insightful. 

Here, we propose a new upper bound using the the derivative of the theta series. 

Theorem 2 (Upper Bound on Distortion). The component Dc is bounded by 

1 / 1 

'z 

where 0^ denotes the derivative of the theta series (with respect to t) 



Dc < -7^0Ae (13) 

27rn V Snai 



Q'Ar) = -Y.''\We-^^\'^\\ (14) 
leA 

Proof: The idea is to upper-bound the integral over Vc by that over the half plane whose 
boundary is equally far from lattice points and 1^, similar to the technique widely used in the 
union bound on the decoding error probability. That is, P(z G Vc(lc)) < > ||lc||/2) = 

Q ^^^^j, where Q{x) = ^^P x > is the Gauss Q-function. Then, we 

have 

Dc < -Y: iiicirg^"'^ 



2n 

IcSAc 
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TABLE I 

Theta series and derivatives of well-known lattices. 





Theta series 


Derivative 








Dn 












Es 


1(^1+^1+^8) 





where the second inequality follows from the fact Q{x) < |exp(— for x > 0. 

Now we may recognize the right-hand side of (fTSi) as the derivative of the theta series. 
Accordingly, the distortion can be expressed in the desired form. ■ 

Remark 1. The bound may be slightly improved by applying the alternative expression of the 
Q-function Q{x) = ^ J^^^ exp (^-^^^ dO for x > 0. 

Now the problem reduces to that of finding a coarse lattice which maximizes the derivative 
of the theta series (note that its derivative is negative). We formalize this new criterion as the 
following Proposition: 

Proposition 2 (Design Criterion of the Coarse Lattice). A coarse lattice is good for Wyner-Ziv 
coding if the derivative of its theta series is large. 

Example 1. Table^shows the theta series and associated derivatives of some well-known lattices. 



C. Even Unimodular Lattices 



As the theta series of even unimodular lattices are polynomial in £"4 and A, it is enough to find 
the derivatives of these two series. The Ramanujan system lfT6l gives the answer. It expresses 
the derivatives of E2, E4 and Eq as functions of themselves: 



dE2{q) 
" dq 

dEi{q) 
" dq 

dEe,{q) 
^ dq 



h (Eliq) - EM) 
\{EMEA{q)-E,{q)) • 
\{E,{q)E,{q)-El{q)) 
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As we are interested in the functions (e ^'^'^) for A; = 4, 6, we get 



y {E2 (e-^-^) E, {e'^^^) - E, (e-^^^) (16) 

^^^^^ = -n {E, {e''-^) E, {e~'-^) - El {e-'-^)) (17) 

and so, combining (fT6l) and (fTTI) . we obtain the derivative of the modular discriminant with 
respect to r: 

= -2nE, (e--) A (e^^) = -^i^. (e^-) {E^ (e^-) - E,^ (e--)) . 

(18) 

This means that, using (|7]), (fT6l) and (fTSi) . we are able to get the derivative of the theta series of 
an even unimodular lattice as a function of E2, E4 and Eq. 

Example 2. The theta series of the Gosset lattice Eg is 

eE,{r)=E,{e-'--). 

Using n6\) . we get 



Example 3. The theta series of the Leech lattice A24 is 

QmM = El - 720 A (e-^^^ = ^ (7^1 (^"^^1 + Si^e (e'^'^l) • 

From this expression, and using ( [7^ and ( (771) . we deduce the derivative with respect to y: 

1) The average behavior: By the Siegel-Weyl formula, the average theta series of an even 
unimodular lattice of dimension 4fc is E2k (e"^'^'^) [[T4l|. This series is a polynomial in E4 and -Ee^ 
and so, its derivative is a polynomial in E2, E4 and E^. We give, here, a method of calculating 
this polynomial. 

Set 

G2k = 2C i2k)E2k 
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and 



4 = (2/c + 3) k\G: 



2k+A- 



Then, dk satisfies the recurrence relation, 




2n + 9 
3n + 6 



(19) 



with do = 3G4 = Q({A)E4 and di = 5Ge = 10C(6)-E'6. From (fT9l) . we can compute all Eisenstein 
series as functions of and Eq. Then, use (fT6l) and (fTTI ) to calculate 'i^^k{e~^'"^)/dT. 

Suppose (t| = 0.01. Figure [2] shows the upper bound on the overall distortion Dn as a 
function of R for some lattices (n = 8, 24, 80), and the average over all even unimodular lattices 
of dimension n = 2400 computed by using the derivative with respect to r of the Eisenstein 
series -E1200 (e"^'"'^). For the 'source coding' component Ds, we use Zador's upper bound on 
Gn [fT3l p. 58]. It can be seen that the upper bound on Dn approaches the Wyner-Ziv limit as 
n increases, meaning the bound is asymptotically tight. 



In this Section, we present a variety of methods to obtain nested lattices and the associated 
simulation results. In our simulation, both X and Y are Gaussian. Also, X = Y + Z with 
Y ~ A^(0, 1) and Z ~ A^(0, 0.01). The nested lattices in our design may or may not be similar. 

A. Clean Similar Sublattices 

Recall that a clean similar sublattice means the boundary of its Voronoi cell does not touch 
any fine lattice points lO. It is known that clean similar sublattices perform better than non-clean 
ones [7J. In this subsection we aim to find clean similar sublattices. A similar technique was 
used in solving multiple-description problem ifTTl . In a multiple-description framework proposed 
in [fTSl . a number of clean similar sublattice constructions were given by extending the result 
in [[8l. We will form our own constructions to be suitably used in nested lattices based on their 
constructions. 

In the two-dimensional case, we introduce a complex-valued multiplying factor ^ = a + bi, 
a,b E Z, which is multiplied to the fine lattice points to obtain the coarse lattice points [[TSl . 
For various ^, this calculation should generate both clean and non-clean similar sublattices. The 



IV. Code Design 
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clean similar sublattices were addressed in IfTSl . where the authors proved the sublattice ^T? is 
clean if and only if the nesting ratio c? + is odd, and the sublattice ^^42 is clean if and only 
if a and h are relatively prime. 

Next, we give the nested lattice construction based on iJ^. We use the ring of Lipschitz integer 
quaternions Hq = {a + 62 + cj + dk^ a, b,c,d E Z}, where i,j,k are unit quaternions. Let an 
arbitrary point in a four-dimensional fine lattice Ap he denoted as Ip = x + yi + zj + wk. Then 
set the multiplying factor ^s ^ = a + bi + cj + dk. Multiplying ^ with the fine lattice points Ip 
gives the coarse lattice points Ic-. The calculated is equivalent to the product shown below: 

^ a —b —c —d 
b a —d c 
c d a —b 
\ d —c b a 

So if the generator matrix of the fine lattice is the identity matrix, the first matrix above can be 
seen as the generator matrix of the coarse lattice. 

More generally, it was shown in [18J that for lattices A = Z"^', / = 1, 2, 3 . . . , if there exists 
a geometrically similar sublattice of nesting ratio A^, then N should be of the form m"/^ for 
some integer m, where n = 4/ is the dimension. From ifTSl . the nesting ratio must be odd 
to make sure the similar sublattices are clean. From our simulation, the nesting ratios which 
are odd perfect squares indeed give better performance than others. The simulation results are 
shown in Fig. |3l 

B. Obtaining Similar Sublattices by Scaling and Rotation 

In this subsection, we present a more general method to obtain sublattices by scaling and 
rotation of the fine lattice. The obtained sublattices are not necessarily clean. We set a parameter, 
(3, as the expanding factor. The way to generate the similar coarse lattice is as follows: 

• First expand the fine lattice by multiplying the expanding factor (3, such that the lattice 
points with the smallest norm greater than zero in the fine lattice become the outermost 
points in the basic Voronoi cell of the coarse lattice. When (3 is an integer, there is always 
a less (or equal) number of coarse lattice points on each shell (formed by points of equal 
norm) of the coarse lattice than in the fine lattice. This requirement can be fulfilled by using 
integer expanding factors. 



/ X \ 

y 

z 

\w / 



(20) 
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• For non-integer expanding factors, some expanded lattice points may not be a subset of the 
fine lattice. This case does not only happen to only one point, but also to a group of points 
with the same norm. The angle between these points and their fine lattice counterparts can 
be detected and calculated. Due to symmetry, these angles for points with the same norm 
are the same. So we can rotate the expanded lattice points by this angle, to match the 
positions of those fine lattice points. Also, the rotation for all the norms should be the same 
since the expanded lattice is similar to the fine one; so the rotated coarse lattice remains 
the similar shape as the original fine lattice. 

The values of (3 giving rise to sublattices are those norms of the fine lattices. So (3 is not 
necessarily an integer. By exploiting tables of theta series in lfT3l . we can find a sublattice. 

Example 4. We give an example in the two-dimensional case for both integer and non-integer 
expanding factors. Let the fine lattice be a hexagonal lattice generated by the following generator 
matrix: 



Then the fine lattice is expanded (and rotated when necessary) to form the coarse lattice 

Mc = P'M.p. For fine and coarse lattices with integer expanding factors (e.g. /3 = 3 and 

(3 = 4:), rotation is not needed. When (3 = \/7, rotation is needed to make the coarse lattice 
points form a subset of the fine lattice, as shown in Fig. E] 

Using this method, we now give the simulation results of the rate-distortion performance for 
nested lattice Wyner-Ziv coding for and A24. 

In the three-dimensional case, we implement the scheme using both D3 and D^. The rate- 
distortion performance is shown in Fig. [51 Compared to the one- and two-dimensional case, the 
three-dimensional scheme gives less distortion. Also notice NLC using Z^g is better than the one 
using D-i. 

The simulation results using D4 and are similar, and are both better than the three- 
dimensional case. At rate 1.5 bits per sample, their distortion gaps to the Wyner-Ziv limit are 
2.52 dB and 1.07 dB, respectively. 

Moreover, for A24, our scheme gives distortion performance closer to the Wyner-Ziv limit 
than the SWC-NQ scheme proposed in [7], especially at low rates. Results in |[7| show a gap 




(21) 
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of distortion 1.53 dB from the Wyner-Ziv limit while our gap is less than 1.28 dB at rate less 
than three bits per sample for A24. See Fig. [6l 

One can also see that the "increasing gap" between the rate-distortion curve of NLC and the 
Wyner-Ziv limit [7] indeed exists in concrete implementation. Nonetheless, from Fig. [6l this 
widening gap as rate increases can be handled by increasing the dimension. In the meantime, 
the rate can not be too high in sensor network applications. So the gap is acceptable with 8- or 
24-dimensional lattices. By increasing the dimension of NLC as the rate increases, a constant 
gap from the Wyner-Ziv limit can be maintained. 

C. Ensemble of Random Nested Lattices 

Code designs above are restricted to those proposed in [fT3l . Although they have acceptable 
performance and low complexity, it is better to have a scheme existing in any dimensions. Hence, 
we use the ensemble of good nested lattice codes proposed in based on the concept of random 
lattices. The random lattice ensemble in ^ can be generated as follows. 

• Take p to be prime. 

• Define a kxn generator matrix M, where Mij is uniformly distributed on (0, . . . ,p — 1), i = 
l,...,fc;j = l,...,ra. 

• Apply Construction A in BU to obtain the lattice A'^. 

The n-dimensional cubic lattice Z" can be viewed as a sublattice of the random lattice A'^ [9]. 
Hence, for any dimension, various nested lattices can be obtained by simply applying different 
linear transformations G to both Z" and A'p. Obviously, the resultant sublattice is not necessarily 
similar. The nest ratio N = p"-~^. 

To make the nested ensemble good for the Wyner-Ziv problem, the fine lattice should be 
good for source coding and the coarse lattice good for channel coding. This is proved in [[T9l by 
extending the results in BH. Also, both the coarse lattice GZ" and the fine lattice GA'^ should be 
good for quantization. To make GZ" good for quantization, G should be the generator matrix of 
a good quantizing lattice [[T3l . e.g., the hexagonal lattice in dimension two and in dimension 
eight. 

However, such lattices are not easy to decode. We use the sphere decoding algorithm described 
in [|20l for the quantization to the random lattices, whose speed is tolerable if n < 24. 
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Examples for dimensions 2, 4, 8 and 24 are given in Fig. |7] The rate-distortion performance 
approaches the Wyner-Ziv limit as dimension increases and the gap is only 1.87 dB within the 
limit for dimension twenty-four. The performance is very close to explicit lattices like the Leech 
lattice. 

Since the elements of M are uniformly distributed over Zj, for large p and since the decoding 
is not easy, this ensemble is of theoretic value only. 

V. Conclusions 

In this paper, we have investigated the rate-distortion function and practical code design 
of Wyner-Ziv coding based on multi-dimensional nested lattices. Under the high-resolution 
assumption, an accurate calculation was developed, and an upper bound expressed in terms of 
the derivative of the theta series was derived. These results can be used to judge the performance 
and serve as a practical guide for choosing good lattices for Wyner-Ziv coding. Several practical 
code designs were presented using multi-dimensional NLC. High-dimensional schemes gave 
performance close to the Wyner-Ziv limit. The performance is even better than the SWC-NQ 
scheme [7J. Compared to [7J, our scheme does not require extra Slepian-Wolf coding based 
on powerful error correction codes, thereby enjoying lower complexity. Hence, the proposed 
scheme may be attractive for applications in sensor networks where simple coding schemes are 
needed. 

This work leaves some open problems. Firstly, there may be room to improve the upper 
bound. Secondly, the derivative of the theta series arising from the upper bound, which is to be 
maximized, is a new problem for lattice researchers. Last but not the least, a more systematic 
approach to low-complexity code design is to be pursued. 
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Fig. 1. Comparison of proposed accurate calculation and the lower bound in I?) for fine lattices and E^. Y ~ N{0, 1) and 
Z ~ iV(0,0.01). 
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Fig. 2. Upper bound on distortion 7?„ versus rate for some lattices (n = 8, 24, 80), and an average lattice (n = 2400). 
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Fig. 3. Distortion versus rate for clean similar sublattices of fine lattices and Z*. 
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Fig. 4. Example of two-dimensional coarse lattice: /3 = \/7. □ coarse lattice points, • fine lattice points. 




Fig. 5. 



Distortion versus rate for fine lattices D3 and -D3. 
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Fig. 6. Distortion versus rate for fine lattice A24. 




Fig. 7. Distortion versus rate for the ensemble of random nested lattices. 



